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NICHOLAS KATZOURAKIS 



Abstract. Let H g C^(K^^") be a Hamiltonian function. Herein we study 
Aronsson Maps, that is solutions u : Q C R" — > to the system 

(1) AacU := [Hp®Hp + H[Hp]^HppyDu):D'^u = 

with emphasis on the special case of 2-d oo-Harmonic maps for n = 2 < N with 
H(P) = ||P|2 the Euclidean norm. We also consider the 1-d case of Aronsson 
ODEs when H {x , u{x) , u' (x)) depends on all arguments. (1) was first derived 
in the author's recent work [I- ']. By establishing a general Rigidity Theorem 
for Rank-One Maps of independent interest, we analyse the phase separation 
of 2-d oo-Harmonic maps and their interfaces whereon the coefHcients of the 
system become discontinuous. As a corollary, we extend the Aronsson- Evans- 
Yu theorem on the non-existence of zeros of \Du\ for solutions to all A'^ > 2 and 
establish a Maximum Principle (Convex Hull Property) for N = 2. We further 
classify all H for which (1) is elliptic: they are the "geometric" ones, which 
depend on Du via the Riemannian metric DvJDu. We also study existence, 
uniqueness and regularity of the initial value problem for the Aronsson ODEs. 



1. Introduction 

Let 7J be a Hamiltonian in C^(M^^"). In this paper we primarily study Aronsson 
maps, that is solutions u : C M" — > to the system 

(1.1) AooU (^Hp®Hp + H[Hp]^Hpp^{Du):D^u = 0. 

In (1.1) [iJp(P)]-L denotes the projection on the nuUspace of Hp{Py : — ^ M" 
and Hp is the derivative matrix (for details see Preliminaries 1.1). Particular 
emphasis is given on the important special case of 2-dimensional oo-Harmonic maps 
ioTn^2< N and H{P) = ^\P\'^, where \P\ = (tr(P^P))i/2 is the Euclidean norm 
on E^^". In this case (1.1) simplifies to 

(1.2) Aoou := (^Du(E) Du+\Du\'^[Du]-^^ : D'^u = 0. 

We also consider the 1-dimensional Aronsson ODE system when iJ (x, u(x), it'(x)) 
depends on all arguments. (1.1) is a quasilinear system in non-divergence form 
which arises in the limit of the Euler-Lagrange equations of the functionals 
H{Du)P as p — >■ cx) and was first derived in the author's recent work [K 'i]. Let 
us briefly recall the derivation in the simpler case of (1.2). After expansion and 
normalization of the p-Laplace system ApU = Div(\Du\P~^ Du) = 0, we have 

\Du\^ 

(1.3) Du®Du:D'^u+^- ^Aii 0. 
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Let [IJw]^ and [Du]-^ denote the orthogonal projections on the range of Du and 
the nuhspace of Du^ respectively. Since [Z?m]^ + [Du]-^ — I, by expanding Au 
with respect to these orthogonal projections, we get 

(1.4) Du(E)Du: D^u + ^- UDu]'^ Au ^ -J UDu]^Au. 

p — 2 P ~ 2 

By orthogonality, right and left hand side of (1.4) are normal to each other. Hence, 
they both vanish and (1.4) actually decouples to 2 systems. By renormalizing the 
right hand side of (1.4) and rearranging, we get 

(1.5) Du<»Du: D'^u+\Du\'^\Du]^Au = ^[Dul^Au. 

p-2 

As p — > oo, (1.5) formally leads to (1.2). In the scalar case of = 1, the normal 
coefhcient |Z?up[Z)u]-'" vanishes identically and the same holds when u is submersion 
map. The scalar oo-Laplacian DiuDjuDjjU = for = 1 has a long history. It 
was first derived and studied by Aronsson in the '60s in [A3, A4] and has been 
extensively studied ever since (see for example Crandall [ ], Barron, Evans, Jensen 
[BE.I] and references therein). A major difficulty in its study is its degeneracy and 
the emergence of singular solutions (see e.g. [A6, A7, Kl]). In the last 20 years the 
scalar PDE has been studied in the context of Viscosity Solutions. 

In the vector case even more intriguing phenomena occur, highlighted in ['••''■]. 
Firstly, there emerge highly singular solutions to (1.2). For example, for any unit 
speed curve / G C^(M)^, the planar map u{x,y) := f{x) + if{y) is a singular 
"oo-Harmonic" local diffeomorphism on a neighborhood of the diagonal. Secondly, 
a further difficulty not present in the scalar case is that (1.1) has discontinuous 
coefficients even for C°° solutions: such an example is given by 

(1.6) u{x,y) := e" -e'y , u : {|x ± ?/| < tt} C _^ 

Indeed, in [i\ we showed that (1.6) is oo-Harmonic on the rhombus, but the rank of 
Du jumps to one on {x — y}, while rk(£'u) = 2 otherwise and hence the projection 
is discontinuous. In general, oo-Harmonic maps present a phase separation. 
On each phase the dimension of the tangent space is constant and these phases 
are separated by interfaces whereon [Du]^ becomes discontinuous. On a phase, we 
interpret (1.2) as decoupling to the tangential system DuD(^^\Du\^) = in 
and the normal system |Z?wP[Z?m]^Am = in [ZJu]^. Much more intricate examples 
of smooth oo-Harmonic maps whose interfaces have triple junctions and corners are 
constructed in [ ]. For further extensions to the subelliptic setting see [K7]. The 
appearance of discontinuous coefficients is a genuinely vectorial phenomenon which 
appears when Tk{Du) > 2. In particular, when n = 1 all oo-Harmonic curves are 
affine and for u : f7 C M — > M^, (1.2) reduces to 

(1.7) A^u = iu'®u')u" + \uf(l-^®^)u" = \u'\'u". 

V \u'\ \u'\/ 

System (1.1) arises as a sort of Euler-Lagrange PDE related to L°° variational 
problems for the supremal functional 

(1.8) Eoc{u,n) -.^ ess sup H{Du). 

n 

In the very recent paper [K4] we identified the variational principle govering oo- 
Harmonic maps for the model functional ||£'u|jioo(Q) of vector- valued Calculus of 
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Variations in Recently, special cases of (1.1) and (1.8) in the vector case have 
attracted substantial interest. Ou, Troutman and Wilhelm in [OTW] and Wang 
and Ou in [\\ ( )] studied Riemannian variants of tangentially oo-Harmonic maps 
which solve only the tangential part of (1.2). Sheffield and Smart in [ ''-] used as 
Hamiltonian the nonsmooth operator norm on 



pNxn 



(1.9) 



IPII 



max (P P : w w) 



and derived a very singular variant of (1.2) for a norm different from the Euclidean, 
which governs optimal Lipschitz extensions of maps. The choice of (1.9) owes to 
that they employ an L°° variational approach and they need the coincidence 



(1.10) 



ess sup 
n 



\Du\\ 



Up{u,n) , ueC\W)^, n convex 



for the Euclidean norm \Du\ of Du on 
Hamiltonian 



pNxn 



Capogna and Raich in [C'R] used the 



(1.11) 



HiP) 



det(P) 



defined for local diffcomorphisms u : C M" — y M" and developed an L°° vari- 
ational approach to optimise Quasiconformal maps. They derived and studied a 
special important case of (1.1). Their results have very recently been advanced by 
the author in [lv.5]. 

Herein we are mostly interested in understanding the structure of system (1.1) 
and of its regular solutions, with emphasis on those of (1.2). To this end, we first 
develop a Geometric Analysis tool of independent interest, and then employ it to 
our analysis. Section 2 is devoted to this analytical result. 

To begin with, consider a map u : il C- M" — > given as composition of a 
scalar function / e C^(il) with a unit speed curve : M — > M^, that is u = i/ o f. 



^(f(0)) 




Figure 1. 



Then, Du = (z^o /) (g) Df and hence m is a Rank-One map, that is tIl(Du) < 1 on f2 
and Du can be written as Du = ^ (X) w for two vector fields ^ and w. Interestingly, 
the class of Rank-One maps is rigid since a certain converse is true as well: all 
maps which satisfy t^(Du) < 1 arise as compositions of unit speed curves with 
scalar functions. More precisely, in Theorem 2.1 of Section 2 we prove that for any 
Rank-One map u £ C^(ri)^ over a contractible domain J7 C M^, there exists an 
/ e C^(17) and a partition of to connected Borel sets Bi such that on each Bi, u 
can be represented as composition of / with a twice differentiable unit speed curve: 
u = o f . Theorem 2.1 is optimal. Without extra assumptions, there may not 
exist any / globally defined on f2 and w(f2) may bifurcate without being given by 
a single-valued curve v for which u — v o f (Corollary 2.2, Example 2.3). 
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Theorem 2.1 is a result ol rigidity type and has been motivated by the rigidity 
results of Rindler in [Rl, R2]. Actually, we extend a part of his result from con- 
stant rank-one tensors ^ iX) w to variable rank-one S.{x) (g) w{x) tensor fields. When 
compared to the rigidity results known in the literature which relate to Gromov's 
Convex Integration (see e.g. Kirchheim [Ki]), Theorem 2.1 is somewhat surprising 
in that most rigidity phenomena appear for rank greater than 2. 

The idea of the proof of Theorem 2.1 has two main steps. Suppose that Du ~ 
^ (g) w. By using that Curl(_Du) = 0, we first show that rk(£'^) < 1 as well and also 
invoke Poincare's lemma of De Rham Cohomology to represent w by the gradient 
Df of a scalar function. Then, we employ geodesic flows, exponential maps of 
Riemannian Geometry and a curvilinear extension of "De Giorgi-type" arguments 
to show that ^ and / locally have the same level sets and hence ^ can be written 
as ^ = j> o /. Hence, locally we deduce Du = {v o f) ® Df — D{v o /). 

It seems that the natural setting for the validity of the Rigidity Theorem 2.1 is 
that of Lipschitz maps. Indeed, we provide such an extension in Theorem 2.4. Yet, 
this does not follow by a direct approximation argument and substantial complica- 
tions arise. The problem is that unless at least one of ^, w is constant, all standard 
mollification schemes seem to fail when one tries to approximate Lipschitz Rank- 
One maps by smooth and Rank-One maps. If both ^ and w vary, the image of the 
mollification may "fatten" and its Hausdorff dimension increases (Remark 2.5). We 
remedy this problem by imposing a rather reasonable extra assumption on Du. 

Using the analytical machinery developed in Section 2, in Section 3 we move to 
(1.2) and analyse the structure of oo-Harmonic maps. In Theorem 3.1 we establish 
that if M : n C R2 — )• M.^ is an cx)-IIarmonic map in C^(0)^, then u naturally 
separates to (at most) two different phases, a 2-dimensional O2 whereon u has 
the structure of an Eikonal immersion and an 1-dimensional fii whereon u has 
the structure of a scalar cx)-Harmonic function along a fixed direction. The two 
phases are separated by an interface S, whereon u has the structure of a scalar 
Eikonal function. Theorem 3.1 relates directly to the phase separation of tight 
maps (vectorial optimal Lipschitz extensions) observed by Sheffield and Smart in 
[SS]. For their non-smooth Hamiltonian (1.9), the 2-dimensional behavior of tight 
maps is Conformal and the 1-dimensional is that of "fans", that is, behavior of 
essentially scalar oo-Harmonic functions but with variable direction. 

By employing Aronsson's theorem in [-v-t], Theorem 3.1 implies an extension of 
the Aronsson-Evans-Yu theorem ([A4, E, Y]) on the non-existence of zeros of \Du\ 
for solutions to (1), for all N > 2. That is, if u : 57 C — y is an oo-Harmonic 
map, either rk(Z?M) > or Du = 0. In view of our earlier observations, solutions 
of (1.2) have positive rank but generally non-constant rank. As a corollary, for 
N — 2 we establish a vectorial version of the Maximum Principle for solutions 
It : C — > M^, known as Convex Hull Property (Corollary 3.3), which states 
that the range is contained in the closed convex hull of the boundary values: 

(1.12) u{n) C c6{u{dfl)). 

The inclusion (1.12) is nothing but an elegant restatement of the Maximum Prin- 
ciple for all projections of u. It is well known in the context of Minimal Surfaces 
(see Colding-Minicozzi II [CM], Osserman [()]) and more generally in Calculus of 
Variations (see the author's paper [K2] and references therein). 

Motivated by Aronsson's paper [AG] and also from ['^'-], in Proposition 3.5 we 
investigate solutions to (1.2) of the radial form u = p'^ f{k9) for fc > and / a curve 
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in M^. Interestingly, such solutions are very rigid, since they always have affine 
range and they moreover are essentially scalar if A: 7^ 1. 

In Section 4 we focus on the general Aronsson PDE system (1.1). We motivate 
our results by observing that (1.2) is quasilinear and degenerate elliptic, that is, for 

(1.13) A{P) -.^ P(g)P+\P\^[P]^(g)I 
we can rewrite the oo-Laplacian (1.2) as 

(1.14) A{Du) : D^u = 

and A satisfies the Legendre-Hadamard ellipticity condition (4.4) and also the sym- 
metry condition (4.5). However, (1.1) is not degenerate elliptic, not even when the 
Hamiltonian H is strictly convex on M^^". The problem is that in the normal part 

(1.15) H{Du)[Hp{Du)]^Hpp{Du) : D^u = 

both tensors [Hp{Du)]-^ and Hpp{Du) are symmetric, but if > 2 their product 
may not commute. Which are the Hamiltonians which lead to elliptic Aronsson 
systems? This question has empty content for N = 1, since Aronsson's equation 
Hp.{Du)Hp.{Du)DfjU = is always degenerate elliptic. In Theorem 4.1 we give 
a complete answer, by classifying the Hamiltonians which lead to elliptic Aronsson 
systems. Every "geometric" Hamiltonian which depends on Du via the induced 
Riemannian metric DvJDu on the range u(Vt) C M^, that is when 

(1.16) H{P) = hilp'^p) , h = h{p), 



.2 

gives rise to a degenerate elliptic PDE system which takes the form 
(1.17) AooM = {ouhp^Duhp + h[Du]^<^hp^ : D^u = 0. 

with h = hiJ^Du^ Duj. In the case of Ago, we have h{p) = tr(p) — pn ... + Pnn- 
In dimensions n < 3, the converse is true as well and this is the only way that 
degenerate elliptic Aronsson systems can arise. However, if n > 4 complicated 
structures in the minors of the higher order derivative tensors -ffp...p(0) appear and 
a necessary extra assumption is required for the converse to be true. Without it, H 
can be written in the form (1.16) up to an 0(|P|'*) correction. In the case n — 1, 
then (1.16) requires that H{P) is radially symmetric: H{P) — /i(i|Pp). This is 
very restrictive, but should be compared with the rigidity of Lipschitz extensions 
for maps in Kirszbraun's theorem (see e.g. Federer [F], p. 201), in contrast to the 
fiexibility of scalar Lipschitz extensions. 

In the light of our general Theorem 4.1, it is not a coincidence that all Hamil- 
tonian known in the literature have the form (1.16). The Euclidean norm trivially 
gives H{P) = itr(P^P). The Hamiltonian (1.9) of Sheffield and Smart evidently 
has this form and the Hamiltonian (1.11) of Capogna and Raich can be written as 



(1.18) H{P) 



/tr(P^P)" Y 



V det (PTP) J ■ 

Moreover, by comparing the elliptic version (1.17) of (1.1) with (1.2), it follows 
that by studying oo-Harmonic maps we gain a great deal of intuition regarding the 
behavior of general Aronsson maps. 

Finally, in Section 5 we focus on the 1-dimensional case of Aronsson ODEs. We 
first formally derive Aronsson's system in the limit as p — >■ c» of the Euler-Lagrange 
equations of an iy^'-functional with Hamiltonian He C^(M x x R^) depending 
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on all arguments: H — H[x,u{x),u' (x)) (equation (5.8)). Then, by arguing as 
in Theorem 4.1, we derive the degenerate elliptic system of Aronsson ODEs with 
dependence on all the arguments: 

(1.19) A^u = \u'\'^(^hpu" - Tlu'hrj^ + h^u' = 0. 

Here h = h(^x,u{x), and R„' is the reflection operator with respect to the 

normal hyperplane bundle [u']-^. As in the case of (1.7), the coefficients of (1.19) are 
continuous, although R„' becomes discontinuous at critical points on {u' — 0}. In 
Theorem 5.2 we study existence, uniqueness and Wi^^{R)^ regularity of solutions 
to the initial value problem for (1.19). 

1.1. Preliminaries. Throughout this paper we reserve n,N G N for the dimen- 
sions of Euclidean spaces and §^^^ denotes the unit sphere of K''^. Greek indices 
a, /3, 7, ... run from 1 to iV and Latin i, j, fc, ... form 1 to n. The summation conven- 
tion will always be employed in repeated indices in a product. Vectors are always 
viewed as columns and we differentiatiate along lines. Hence, for a,b G E", a^b is 
their inner product and ab^ equals a ® b. If u = UaCa '■ ^ C M" — y M.^ is a map, 
the gradient matrix Du is viewed as DiUaCa^ei : — > M^^" = M^(g)E" and the 
Hessian tensor D'^u as Dl^UaCa^ei^ej : fi — > ]R^(g)S(IR"). If F is a vector space, 
then E>{V) denotes the symmetric linear operators T : V — > V for which T = . 
The Euclidean (Frobenious) norm on M^^" is |P| = (PatPai)^ = (tr(P^P))5. If 
F g C"^(IR^^") is a function and we denote the standard basis elements of M^^" 
by Cai := Ca (i) Ci, then its q-ih order derivative tensor Fp p at Pq 

(1.20) Pp...p(Po) =i^p„,.,...p.,.,(^'o)eain ®-®e„,,, 

is viewed as a multilinear map — > (8)(«)E^, or equivalenly as an element 
of (g)(«)(M" ® M"). Here "(g)(«'" is the g-fold tensor product. Hence, Pp...p is a map 

jjATxn ^ 0(9)(MJVxn)_ ^j^^^^ ^ ^^^^^ ^-g^g^j. (j ^ (^i-i) {R^ xn^^ jg 

fully symmetric in all its arguments when 

(1-21) C ...ai...fij... = C,,,aj...j3i... = C , , ,pj ,ai... ■ 

We also introduce the following contraction operation for tensors which extends the 
trace inner product P : Q = tr(P^Q) = P^Qai of M^""". For, if C € ®(«)(IR^><") 
is a g-th order tensor and A € (g)(p)(]R-'^^") is a p-th order tensor with p < q, we 
define 

(1.22) C:A := (C„,,,...„,,,^„^,^...„,. Je„,,, ® ... ^ e„^+,,^^, € (M^>^"). 

Let now P : E" — ^ be linear map. Upon identifying linear subspaces with 
projections on them, we have the split = [P]^ © [P]^ where [P]^ and [P]^ 
denote range of P and nuUspace of P^ respectively. Hence, if C G §^^^, then is 
(the projection on) the normal hyperplane /— Consequently, the cx)-Laplacian 
(1.2) in index form reads 

(1.23) D.Uo^DjUpDlup + \Du\^[Du]ipDlup = 
and the Aronsson system (1.1) becomes 

(1-24) {Hp^.Hp^^ + H[Hp]i^Hp^^p^^){Du)Dlup = 0. 

For convenience we use a different scaling in (1.23) and (1.24) and we multiply the 
normal term of (1-23) by a factor 2 which is plausible since (1-24) consists of two 
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systems normal to each other. Finahy, 'h'^ denotes the fc-dimensional Hausdorff 
measure and for measure theoretic notions we use herein we refer to Simon [S]. 

2. Rigidity of Rank-One Maps. 

2.1. The case of smooth Rank-One Maps. In this subsection we establish our 
Geometric Analysis rigidity result in the case of maps. 

Theorem 2.1 (Rigidity of Rank-One Maps). Suppose ft C M" is open and con- 
tractible and u : fl C_ M" — y M.^ is in C^(fi)^. Then, the following are equivalent: 

(i) u is a Rank- One map, that is rk{Du) < I on or equivalently there exist maps 
^■.n — >R^ andw.n — ^ with w G Ci(r2)" and ^ e C^{n\{w = 0})^ such 
that Du = ^ (g) w. 

(ii) There exists f G C'^{fl), a partition of Q. to Borel sets where each Bi 
equals a connected open set with a boundary portion and Lipschitz curves {i^'jigN ^ 
Wi^^ (M.)^ such that on each Bi u equals the composition of v"^ with f: 

(2.1) u ^ v'o f , onB.Cfl. 

Moreover, |z>'| = 1 on f(Bi), z>* = on M \ f{Bi) and there exists i)^ on f{Bi), 
interpreted as 1-sided on df{Bi), if any. Also, 

(2.2) Du ^ {v^ o f)®Df , on B, C Vl, 
and the image u{^) is an 1-rectifiable subset ofM.^: 



(2.3) u{n) = U^''(/(^0) c 



As we have already mentioned in the Introduction, an extra assumptions is 
required in order to deduce that a rank-one map u has the form u = vof for a unique 
single- valued unit speed curve v. This assumption guarrantees "low complexity" 
for the direction field f . 

Corollary 2.2 (Strong Rigidity of Rank-One Maps). Suppose Q, C M" is open 
and contractible and u : C_ M" — y M.^ is in C^{il)^ . Consider the following 
statements: 

(i) u is a strictly Rank- One map, that is rk{Du) = 1 on or equivalently there 
exist maps ^ : O — > \ {0} and w: n — > R" \ {0} such that Dm = ^ (g) u>. 
Moreover, the following condition holds 



(2.4) E 



nn (^[j d{\DU > 0}^ = 



(a) u equals the composition of a single curve v g {R^^ with a scalar function 

f G C'^{fl), without critical points that is u = v o f with \v\ = 1 on f{^), i> = Q on 
R \ fi^t). Moreover, Du ~ (i' o /) ® Df on f2 and u(f2) is l-rectifiable, equal to 

Then, (i) implies (ii) and also (ii) implies that u is a strictly rank-one map, that 
is assertion (i) without (2.4). 
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Example 2.3. The additional assumption (2.4) of Corollary 2.2 is necessary in 
order to obtain u — v o f . It reduces the complexity of ^ and leads to the avoidance 
of bifurcations in the curve v. For, let u :I 

{ + fHx)Jix)y, 

u{x) := <j f\x)J\x)Y , 



— > ]R~ be given by 
on {/ > 0}n {a;i > 0}, 
on {/>0}n{a;i < 0}, 
{/ < 0}, 



where 



fix) 1 - \x-ei\^\x + ei\' 



Then, u can not be written as u ~ v o f for a single-valued curve v since the unique 
V bifurcates and has two branches: i^^{t) = ( ± <^X(o,oo)(^)7 ^) • 




Figure 2. 



Proof of Theorem 2.1. The imphcation (m) ^ (i) is trivial and the whole proof 
is devoted to establish the reverse implication (i) (ii). For, suppose there exist 
^ : C — s- and w : C R" — s- M" such that Du = i®w. By replacing ^ 
by C/ICI on {|^| > 0} and w by \£^\w on {|^| > 0}, we may pass all the zeros of Du 
to w and assume that |^| = 1 on 

(2.5) Q.0 {|Du|>0} = {\w\>Q}. 
By differentiating DkUa — S.aWk, we have 

(2.6) D^u^ = iD,Uw^ + UDjW,), 



(2.7) 



Since u G C^(ri)^, the curl of Du vanishes and we have 

(2.8) D^ju^ = D%u^. 
Hence, by (2.6), (2.7), (2.8), 

(2.9) iD,Ow^ - {D,i)w, = i{D,Wj - Djw,). 

Since 1 on fioi we have I?fc^^^ = thereon. Hence, the two sides of (2.9) 



are normal to each other. By applying the projections ^ ® ^ and 
(2.9) decouples on JIq to 



(2.10) 
(2.11) 



Curi(u;) 



DiW. 



DjWi 



0, 



0. 



By (2.10), the curi of u; : f^o ^ M" — > vanishes and by (2.5) w = on rj \ f^o ■ 
Hence, since is contractible, by Poincare 's Lemma w can be represented by the 
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gradient of a scalar function / G C^(il): w — Df. By (2.11), for all i,j e {1, n} 
for which {wi ^ 0} C] {wj 7^ 0} 7^ 0, we have 

(2.12) = 

Wj Wi 

By (2.12), the quotient -Dfc^a/wfc is independent of k. Hence, we may define 

D ^ 

(2.13) 77 — : {wk 7^ 0} C — > M^. 

By (2.12), 77 is well defined on all of ilo since U"{K;fe 7^ 0} is an open cover of 
Hq — {\w\ > 0} and on the overlaps the different expressions coincide. By (2.13), 
we have Dk£,a — VaWk on {wk 7^ 0}. Actually, this extends to the whole of fio since 
by (2.11) we get Dk^ = whenever Wk = 0. Thus, 

(2.14) D( = rj®DJ , on f^o, 

and also 77 is normal to ^, since 77^^ = ^^-DfcdlCP) — on {wk 7^ 0}. We now 
employ (2.14) to show that in a certain local sence ^ and / have the same level sets. 
Fix a e {1, A^} and set 

(2.15) A r!on{|77„| > 0}, 

(2.16) g := ^„ , A := 77^. 

We then obtain 

(2.17) Dg = XDf , on A, 

while \Dg\ > and |A| > on A. (2.17) says that the level hypersurfaces {/ = f{x)} 
and {g = g{x)} passing through x have for all x G A the same tangent spaces: 

(2.18) [Dg]^ = [Df]^ = I ^ 



\Df\ ^ \Df\- 

Consider the level hypersurfaces of /, g as Riemannian submanifolds of A with the 
induced metrics from K". Since covariant derivatives coincide with tangential pro- 
jections of derivatives in K", the geodesic equations for x, "0 with initial conditions 
X(0) = -0(0) = a; e ^ and x(0) = V'(O) = e e [Df{x)]^ = [Dg{x)]^ are 



(2.19) 



[Dfixim^m = 0, t>o, 

X{0)^x, x(0) = e, 



I [Dgim)]H{t) = 0, t>0, 

Since [i'.g]"'" = [Df]-^, x and tp satisfiy the same ODEs with the same initial condi- 
tions. Hence, by uniqueness, % = 7/;. Consequently, the esponential maps exp^ and 
expg of {/ = f{x)} and {g — g{x)} coincide and hence (expg)~^ o exp^ equals the 
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identity their common geodesically convex neighborhood centered at x. 




Figure 3. 

Hence, the level hypersurfaces of /, g within A coincide, but perhaps they are at 
different heights. Cover A by countably many balls whose radii are small enough 
to guarrantee that the intersections of the level sets of /, g with each ball are 
connected. Using this cover, we decompose A to a partition of connected Borel 
sets by writting A = W^Ai, where each Ai equals an open subset of the ball of the 
cover with possibly some boundary portion. Then, for each i G M and each i G N 
there is a unique G K such that {f = t} equals {g = p^{t)} locally within Ai. 
Hence, there exists a unique bijection p* : f{Ai) C M — > giAi) C M such that 

(2.21) {g^P^t)} = = o = {p'of^p\t)}, 

within C fig. Equivalently, 

(2.22) g = pi o f , on A„ i e N. 
We extend from /(Aj) to M by zero. 

g(Ai) 




Figure 4. 



On \ ^ = ^0 \ U^Aj, we have Dg = 0. Hence, there exists a constant function 
p"" ■■ fi^o \A)CR- 



such that 

9 - p"of 



on \ U^A,. 
as well. By recalling (2.15) and (2.16), we have shown 



(2.23) 

We extend p*^ by zero on 
that for any ^a, 1 < < A^, there exists a partition of flo to disjoint connected 
Borel sets Af where each Af equals an open set with possibly some boundary 
portion and also their complement Aq := \ Uf^Af. There also exist functions 
Pa'.R — >R such that 

(2.24) Co = Piof , on A?, z = 0,1,2,... . 

Hence, by recalling that |^| = 1 on IIq, there exists a partitition of Jig to connected 
Borel sets which are intersections of the Ai^s and respective bounded curves 

p* : M — ^ {0} U S^-i C which satisfy 

(2.25) |p*| = 1 on f{B,) , p^ = on M \ f{B,), 
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and are such that 

(2.26) ^ = o / , on B„ 
for alH e N. We set 

(2.27) i^'it) := f n\s)ds , i £ 

Ja 



N. 



Then, by (2.25) we have that e W{-^^(R)^, while |z>*| ee 1 on the interval f{B,) 
and also = on M \ f{B^). By (2.26) we have 

(2.28) ^ = i^'of , on B,. 
Hence, (2.28) implies 

(2.29) Du = ^(^w = {u''of)®Df = D{iy'of), 

on Bi. Thus, u = o f on each Bi C fig, up to an additive constant. By taking 
difference quotients in (2.28), comparing with (2.14) and passing to limits, we obtain 

(2.30) DC = {i)'of)®Df, 

and hence i^a ° f = Dk^a Dkf , on Bi. Thus, i/^ exists on f{Bi) C M and is 
interpreted as 1-sided at the cndpoints of this interval in case it is not open. Since 
Du = and Df = on 9(0o) n f2, we can extend the partition VJf'Bi of ilo to 
f2o ri f2 and further extend the families {-BijigN and {t^*}igN by attaching the limit 
values and setting 

(2.31) So := f7\no, 

(2.32) := = const. 

Hence, since u — i/^ o f on each Bi of the partition Ug°i?i = 51, we conclude that u 
is 1-rectifiable and the image u{Q) equals a union of images of Lipschitz curves: 

CO 

(2.33) u{n) = |Jz.'(/(i?,)). 

i=l 

The theorem follows. □ 

Proof of Corollary 2.2. In the setting of the proof of Theorem 3.1, if in addition 
the set E given by (2.4) is empty and moreover rk(Du) > on fi, then for all a G 
{1, ...,7V}, either DS,a does not vanish anywhere inside JIq = fl or it is identically 
constant. In both cases, the previous set A is connected and coincides with Q. 
Hence, the curve v constructed is unique and consequently u = v o f with |!>| = 1 
on f{ft) and j> = on M \ f{il). The reverse implication is obvious. □ 



2.2. Extension to Lipschitz Rank-One Maps. In this subsection we extend 
Theorem 2.1 to the Lipschitz setting. As we have already explained, this does not 
follow by a straightforward mollification argument and an additional assumption is 
required. 

Theorem 2.4 (Rigidity of Lipschitz Rank-One Maps). Suppose ft C M" is open, 
bounded and contractible and u : fJ C M" — y M.^ is in W^'°° {il)'^ . 
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We moreover assume that there exists a family {V^'^}e>o of rank-one smooth 
tensor fields in C°°(f2)^" where each is curl-free (that is rk{V'^) < 1 and also 
DjV^, - D,V^j = O; such that 

(2.34) V Du in L°°{n)^" and V — > Du a.e. on fl, as e ^ 0. 
Then, the following are equivalent: 

(i) u is a Rank-One map, that is rk(Du) < 1 a.e. on fl or equivalently there exist 
L°° vector fields ^ : VL — > and w : Vl — > M" such that Du = ^ a.e. on ft. 

(ii) There exists f G W^'°^{Q,), a partition {-B^jigN of Q. to measurable sets which 
covers it a.e., that is \Vt \ (U^i3i)| — and Lipschitz curves C W^^^ {M.)^ 
such that on each Bi u equals the composition of with f: 



Remark 2.5. The extra approximation assumption (2.34) of Theorem 2.4 requires 
that Du is in the intersection of the weak* and the pointwise closures in L°°(ri)^" 
of the cone which consists of smooth rank-one curl-free tensor fields. Such an 
assumption is superfuous if either ^ or w is identically constant, since mollification 
of Du = £,^vu produces the desired approximations . 

Generally, however, all standard mollification methods average at each point 
contributions from nearby points. As a result, if such a "partial affinity" of u fails 
to hold and both ^ and w vary, the range u{n) may "fatten" and the mollification 
of u may not be rank-one any more. 

Unfortunately, we have not been able neither to verify the necessity of the as- 
sumption nor to construct a proper mollification scheme allowing to drop it. 

Proof of Theorem 2.4. Is suffices to demonstrate the implication (i) (m). 
Suppose Du = ^ w a.e. on fi. By a rescaling of the form Du ~ ® (ICI""^) 

{1^1 > 0}, we may assume that ^ : Hq — > S^~^, where := {\Du\ > 0} C ft and 
also that £, = a.e. on r2\ Hq. By assumption, we have rk(F'^) < 1 and hence there 
exist : C M" — > and -.ncW — such that ^ (g) . By an 
appropriate rescaling inside the products (j^^^) (|C^|t«'^) on {|^'^| > 0}, we may 
assume that ^ : il^ — > S^'^ where := {{Vl > 0} C n and also that = 
on J7 \ Slg. 

We now claim that ^'^ — > ^ and also that w'^ — > w as e — ?> 0, both weakly* 

in L°°{fl) and also a.e. on il; indeed, there exists rj such that — ^ rj and hence 
by the L^(ri)^" strong convergence of (E> which follows by the Dominated 
Convergence theorem, we have 



(2.35) u = I/' o / , onB.Qn. 
Moreover, |j;>*||ioo(R) < 1 and = a.e. on M \ f{Bi). Also, 

(2.36) Du = {v'of)®Df , a.e. onB.cn, 
and the image u{^l) is an 1-rectifiable subset ofR^: 



(2.37) 




(2.38) 



= (f )^(f 0w^) 
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as £ — > 0. Thus, by uniqueness of limits of £,^^w^ we have [{ri(E)r])£^]^w — a.e. 
on D, and hence ^ = rj. Since Q is contractible, by Poincare's lemma, for any e > 
there exists a smooth map : ^l C M" — y such that can be represented as 
the gradient of u"^: Du'^ = (g) w^. Moreover, each is a smooth rank-one map: 
by Theorem 2.1, there exist scalar functions £ C°°(f2), partitions of fl to Borel 
sets {Bf }»eN with ft = Uj°Sf , families of Lipschitz curves {l''''}^efi C 
with ||!>*^||i=o(]{) < 1 and i>'^ = on M \ f^{Bf) such that = o on each 
Bf C £7, while the images u'^{n) are 1-rectifiable, equal to Uf^i'^^{f^{Bf)). 

We will now show that appropriate normalized shifts of the maps approximate 
u. Fix a point x £ il and set d := diam(ri). Since Du'^ — ^ Du in i°°(r2)^" as 
£ — ?> 0, for all x,y £ il and £ > small we have 

(2.39) \u'{x)-u'iy)\ < {\\Du\\Lo.^n) + l)\x - y\. 

We further normalize by considering appropriate shifts, denoted again by u'^ , 
such that u^(x) — u{x). By (2.39), we have 

(2.40) IkllL-(n) < d{\\Duho.^n) + l) + \u{x)\. 

Hence, there exists v such that u'^ — in W^'°°{ft)^ as £ — >■ 0. We will now show 
that u = V. Since Du'^ — > Du a.e. on fl, for H"^^-a.e. direction e S S"~^, we 
have that Du'^ — > Du H^-a-e. on the set {x + span[e]) n =: /. We fix such an 
e. By Egoroff's theorem, for any a G (0, 1), there is an "H^ -measurable set Ea C / 
with H^{Ecr) < (J such that Du^ — > Du uniformly on / \ as £ — 0. Since 
u'^ix) = u{x), by the 1-dimensional Poincare inequality, for £ > small we have 

d pd 

\u'^ {x + te) — u{x + te)\ dt <d \Du'^ {x + te)e — Du{x + te)e\ dt 



'0 JO 

(2.41) < sup \Du' ~ Du\ + d{2\\Du\\Lo.(<_^) + l)n'^{E„). 

I\E„ 

Since u'^ — > v in C°{n)^ and Du^ — > Du in C°(/ \ E^)'^" as £ ^ 0, by passing 
to the limit in (2.41) we obtain 

(2.42) / \v{x + te)-u{x + te)\dt < d (2||i:»it||i.,o(f2^) + l)cr. 
Jq 

By letting ct — ?> 0, by (2.42) we get u = v on I C ft. Since this holds for H"~^-a.e. 
direction e G §^^^, we get it = w on ft. Hence, — ^ u in W^'°°{ft)^ as £ — > 0. 
Since Df^ — ^ w in L°°(S1)" , for £ > small we have 

(2.43) \.r{x)-riy)\ < (||u;|U»(f,)-f l)|x-2/|. 

We further normalize the family by considering appropriate shifts denoted again 
by such that /'^(a;) ~ f{x)- By replacing also each v" with the translate 
z^^*(_ — {fix) — f^ix))), we do not affect the previous normalization u^{x) — u{x). 
Consequenly, (2.43) implies 

(2.44) < d(lkllL==(a) + l) + l/(5?)|. 

As a result, there exists an / such that — ^ / in W^'°°{f^) as £ — > 0. 
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Since ly" ° — on Bf and z)"^* o = on \ Bf, for e, (5 > small we have 

\B!ABf\ = / Ixs^-XbH 
Jn 

(2.45) = / ||i>^^o/-|-|;>*'o/^|| 

Jo 

Since — >■ ^ in L^(ri)^, for each i G N the family {-Bf}e>o is Cauchy in measure 
and hence has a measurable limit i?j C fj. Since for all e > we have f2 = U^Sf 
and Bf O Bj — % ioi i ^ the limit family forms a cover of except perhaps for 
a nuUset: |0 \ (UfBi)\ = 0. We recall that we have = o on and also 
li'>^lL-(K) < 1 and ;>^' = on M \ ^ (B,f). 











rnBfH 




+1 


I 1 1 ! /,si 




-1 


Figure 5. 



Hence, if x G Sf , for any t e M we have 

(2.46) = it-/(^)i + wm- 

lix^ Bf , then f{x) is in the complement of the interval /'^(_B|) and since is 
constant on M \ f^lBf), for any i S M we have 

(2.47) < |t-/(5;)| + max|i.-(9(/^(i?f))| 

< \t-fix)\ + h1|L~(0). 

As a result, since the family is uniformly bounded on Q, for each i G N the 
family {iy^'}e>Q C 1^/^'^°°(M)^ has a weak* limit i/' € VF/j^°°(M)^ which satisfies 
||^''IU°°(R) — 1- By passing to the limit as e — )■ we get u — v'^ o f on Bi C Q and 
= on M \ f{Bi). Finally, the image u{il.) is 1-rectifiable in and up to an 
-H^-nuUset of M^, we have u{n) = Ufi/*(/(B,)). The theorem follows. □ 

3. The Structure of 2-Dimensional cx)-Harmonic Maps. 

In this section we use the Rigidity Theorem 2.1 of Section 2 to analyse the phase 
separation of classical solutions to (1.2) when n — 2 and N > 2. 

Theorem 3.1 (Structure of oo-Harmonic Maps). Suppose u : il C — y , 
where u is an oo-Harmonic map in C'^{i})^ , that is solution to 

(3.1) Aooti = (^Du(g) Du + \Du\'^[Du]-^(E) : D'^u ^ 0. 

Then, there exists disjoint open sets ^1,^2^^ o,i^d a a closed nowhere dense set 
S such that = Qi U S U Q2 and: 
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2 and the map u\^^ 
const. > 0, 



Eikonal 



is given by an 



(i) On 0,2 we have rk(Du) = 2 and the map : 
immersion with 

(3.2) iDup 
on each connected component of 5^2 • 

(ii) On Qi we have rk{Du) — 1 and the map u\ 
essentially scalar oo-Harmonic function f G C^(rii); 

(3.3) u ^ a + ^f , A^f = 0, aeM^, ^eS""-', 

where ^ and a may vary on different connected components o/fii. 

(Hi) On S, \Du\'^ is constant and also rk{Du) — 1. Moreover if S = dQ,i H dQ,2 
(that is if both the 1— and 2— dimensional phases exists, one on each side of the 



interface S) then . S — 
(3.4) u ^ a + U 



is given by an essentially scalar Eikonal function: 
\Df\ = const. >0, aeM^, ^ e 




Figure 6. 

By employing Aronsson's result in [ ] on the non-existence of zeros for the gradient 
of scalar oo-Harnionic functions on the plane, Theorem 3.1 readily implies the 
following 

Corollary 3.2 (oo-Harmonic Maps have Positive Rank). Let u : QM? — !• 

he an oo-Harmonic map in C'^{n)^ . Then, either \Du\ > on ft or \Du\ = on 
O. Hence, non- constant oo-Harmonic maps have positive rank. 

Proof of Theorem 3.1. We begin by assuming N > 2 and setting 

(3.5) ni := int{rk(Du) < l}, 

(3.6) n2 ■■= {rk(i:>M) = 2}, 

and let also S := il \ {ili U 5I2). Our PDE system (3.1) decouples to 

(3.7) Dud{^\Du\^) = 0, 

(3.8) \Du\^[Du]^Au = 0. 

On fl2, we have rk(_Du) — 2 and hence : il2 — > is an immersion. Thus, 
Du{x) possesses a left inverse {Du{x))~^ for all x E O2. Hence, (3.7) implies 

(3.9) {Du)-^DuD(^^\Du\^^ ^ 

and hence d(^^\Du\^^ = on fl2, or equivalently 

(3.10) \Du\^ = const., 
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on each connected component of $72- Moreover, (3.10) holds on S as well, the 
common boundary of and ili. 

On the other hand, on fii we have ik{Du) < 1. Hence, there exist vector 
fields ^ : f]i C ]R2 — ^ and w : f)i C — ^ ^^^^i that Du = w. 

Suppose first that fii is contractible. Then, by the Rigidity Theorem 2.1, there 
exists a function / g C^(ili), a partition of fli to Borel sets {Bi}i^m and Lipschitz 
curves {iy'}ieN Q Wl^^{R)^ with |z>'| = 1 on f{Bi), = on M \ /(b/) twice 
differentiable on f{Bi), such that u = I'^of on each B^ and hence Du ~ {v'^of)®Df 
on Bi. By (3.7), we obtain 



((;>' o /) ® i?/) ® {{v^ o /) ® Df) 



= 0, 



(3.11) : [(i>'' of)®Df®Df + [v' o /) ® 

on i?i C ill. Since |!>*| = 1 on f{Bi), we have that is normal to and hence 

(3.12) ((;>■' o /) ® Df) ® ((z>^ o /) ® i^/) : ((;>' o /) ® i^^) - 0, 
on Bi C r^i. Hence, by using again that |j>*p = 1 on f{Bi) we get 

(3.13) {Df®Df:D^f){i^'of) ^ 0, 

on i?i C 0,1. Thus, Aoo/ = on Bi. By (3.8) and again since |!>*p = 1 on f{Bi), 
we have [Z?u]-'- = o /]-'- and hence 

(3.14) \Df\^[i''of]^Biv{{i^'of)<E>Df) = 0, 
on Bi Ci^i. Hence, 

(3.15) \Df\' [v^ o f]^({P o f)\Df\^ + (;>» o /)A/) = 0, 

on Bi, which by using once again |;>*p = 1 gives 

(3.16) \Df\\v'of) = 0, 

on Bi. Since Aoo/ = on and f^i — UfBi, f is cx)-Harmonic on fti. Thus, by 
Aronsson's theorem in [A4], either \Df \ > or \Df \ = on fti. 

If the first alternative holds, then by (3.16) we have = on f{Bi) for all i 
and hence is affine on f{Bi), that is i^'(t) = + a' for some |^'| — 1, a' S K^. 
Thus, since u — v"^ o f and u G C^(rii)^, all and all a' coincide and consequently 
u = ^/ + a, C G S^"\ where a G and / S C2(f7i). 

If the second alternative holds, then / is constant on i7i and hence by the 
representation u = ly^ o f, u is piecewise constant on each Bi. Since u G C^(f7i)^ 
and Oi = UfBi, necessarily u is constant on fli. But then IDmIo^I = l-D/ls] = 
and necessarily = 0. Hence, \Du\ = on $7, that is u is affine on each of the 
connected components of fl. 

If r^i is not contractible, cover it with balls {BmjmGN and apply the previous 
argument. Hence, on each B^, we have u = C"*/™ + a™, G S^-\ a™ G and 
/" e C^{Mm) with Aoo/" = on B„ and hence either |D/™| > or |D/™| = 0. 
Since u G C^(J7i)^, on the overlaps of the balls the different expressions of u must 
coincide and hence we obtain u = £,f + a for ^ G S^~^, a G M.^ and / e C^(Oi) 
where ^ and a may vary on different connected components of f2i. The theorem 
follows. □ 

Theorem 3.1 implies a vectorial version of the Maximum Principle when n = 
N — 2, which we now prove. 
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Corollary 3.3 (Convex Hull Property). Suppose that u : C — y is an 
oo-Harmonic map. Then, for all ^l' CC fi, the image u{U,') is contained in the 
closed convex hull of the boundary values: 

(3.17) u{n') C co{u{dn')). 

Proof of Corollary 3.3. We begin by observing that (3.17) is an elegant restate- 
ment of the Maximum Principle for all projections tj^ u of u, that is, when for all 
il' CC and all directions rj e S^~^ we have 

(3.18) sup?7^M < max?7^M. 

Indeed, (3.18) says that w(il') is contained in the intersection of all halfspaces 
containing u{dfl'). To see (3.18), fix ^l' and rj € S^~^ and let ili, fl2, S respectively 
be the constant rank domains and the interface of u, as in Theorem 3.1. Suppose 
that u = ^/ + a on f2i U 5, where C € S^"\ a G and / G C^ini U S). Then, 

\D{tj'^u)\ = \v'^Du\xn, + Iv'^ Du\xsun, 

(3.19) - \v'^Du\xa, + \v'^ ^\\Df\xsun,. 

If \Df\ = on ill, then — ^ and u is affine. Hence, (3.18) follows. Suppose now 
\Df \ > on ill. Since u\fi^ is a local diffeomorphism, we have \ri^ Du\ > for all 




Figure 7. 

Consequently, for all rj E E>^^^ \ [S,]-^, in view of (3.19) we have \D{t]^u)\ > on 
il. Hence, rf^u has no interior critical points inside 51 and consequently we have 

(3.20) niax?7^u = maxjy^u, 

for all directions ?y / 6 By letting dist(?7, [(]-^) 0, (3.20) implies (3.17). □ 

As a topological consequence of the Convex Hull Property for planar cx)-Harmonic 
maps, we deduce the following: 

Corollary 3.4 (Absence of Interfaces). Suppose that m : 51 C — y is an 

oo-Harmonic map. Then: 

(i) If fl2 ~ {rk{Du) = 2} CC Jl, then fl2 = o,nd 5 = 0. Hence, either the set 
whereon u is a local diffeomorphism has a common boundary portion with 51 or it 
is empty and u is everywhere essentially scalar without any interface S . 

(ii) If CC and u extends on dil as an essentially scalar function, there is no 
interface S inside 51 and u is an essentially scalar oo-Harmonic function on 51. 
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3.1. Rigidity of Radial 2-Dimensional Solutions. In this subsection we study 
a class of special solutions of the oo-Laplacian, that of smooth cx)-Harmonic maps 
, TV > 2 of the form u = p^f{k6) in polar coordinates {p,0). Here 



pN 



fc > is a parameter and / : M — > R is a curve in M . It follows that such 
solutions are very rigid, because if fc 7^ 1 they are essentially scalar and if fc = 1 
they always have affine image. Accordingly, we have 

Proposition 3.5 (Rigidity of Radial 2-D oo-Harmonic maps). Suppose that u : 
— > M^, N > 2, in an oo-Harmonic map of the form u = f{kO) in polar 
coordinates (p, 6*) G fc > 0, / G C°°{R)^ . Then, f solves the ODE systems 



(3.21) 
(3.22) 



/'»/'(/" + /) 



1 



(l/f + \f?)f 



Moreover: 

(i) If k 1, then all such solutions have constant rank one and the image u(]R^) 
is contained into a line passing through the origin and f can be represented as 
= i9{d) for some ^ G §^"1 and g G C°°(M). 

If k = 1, then all such solutions have rank at most two and the image u(]R^) is 
contained into a 2-plane ofR^ passing through the origin. On this plane f can be 
represented by 

(3.23) f{e) = ccosB(6l) (cosA(6'),sinA(6'))^, 

where c G M and A, B e C°°(M) satisfy the identity \B'\^ + jA'pcot^B = 1 and 
also < B < f . 

Proof of Proposition 3.5. The derivation of the "tangential part" (3.21) of Aoo 
is entirely analogous to Aronsson's derivation of its scalar counterpart in the paper 
[A6], p. 138. Hence, it suffices to outline the derivation of the "normal part" (3.22). 
Since for all a G {1, ---jN} we have Ua — p^fa{k6), we obtain 



DyUf^ 



(3.24) 



cos ( 
sin^ 

cos ( 
sin^ 



— smi 

cos 9 

— sin ( 
cos 6 



DpUa 

kp'^-'fo. 
kp'-'fL 



Let O{0) denote the rotation-by-0 appearing in (3.24). By recalling that {f,f') is 



pN 



we may write 



a matrix- valued curve M - 

(3.25) Du = kp''-\f,f)0{9)\ 
Hence, since 0{9)^ = 0{9)^^ we have 

(3.26) N{Du'') = {,7GM^ : r/^ (/,/') 0(0)^ - 0} ^ N{ifJ')'^). 
and consequently [Du]-'- = [(/, /')]^- Moreover, 



-DpU 



(3.27) 



[Du]^Au =[(/,/')]^ 

= [ifjT{kp''-'f + k{k-l)p'^-\f + k'p'-'f) 
= fc2/-2[(/,/')]^r. 
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By Corollary 3.2, we may require \Du\ > and hence (3.22) follows by (3.8) and 
(3.27). Now, for (i) we have that if fc =^ 1 then on {|/| > 0} (3.21) gives 

Consequently, /' is everywhere proportional to / and as a result /(M) is contained 
into an l-dimensional subspace of M^. 

For (a), we have that if fc = 1 then (3.22) implies /" = A/ + ///' for some A, /.t G 
C°°(]R). Hence, /(K) is contained into a 2-dimensional subspace of M^. (3.21) gives 
the extra condition that f'^{f" + /) = which implies + |/p = for some 
c G R. Hence, if c ^ then ^(|/'|, |/|)^ is on the unit circle and as such |/| = ccosi? 
and l/'l = csini?, for some B valued in [0, |]. Hence, / = ccosi3(cos A,sinyl)^ 
for some A. The differential relation + cot^ B = 1 follows easily. □ 

4. Classification of Elliptic Aronsson PDE Systems. 

In this section we focus on the general Aronsson system (1.1). As already ex- 
plained in the introduction, when N > 2 the normal coefficient H[Hp]-^Hpp is not 
commutative and as a result the system generally is not degenerate elliptic, not 
even for strictly convex Hamiltonians. In Theorem 4.1 below we establish that all 
"geometric" Hamiltonians which depend on Du via the induced Riemannian metric 
DvJDu lead to elliptic systems. Moreover, in low dimensions n < 3 the converse is 
true as well for (normalized) analytic Hamiltonians with fully symmetric Hessian 
tensor. When n > 4, there appear complicated structures in the minors of forth and 
higher order derivatives and an additional assumption is required. The constructive 
method of proof reveals that it is necessary. 

The main idea in the reverse direction is to impose the commutativity relation 
[Hp]^Hpp = HpplHp]-^ and use power-series expansions of H to derive the form 
(1.16) inductively, by a term-after-term blow-up argument along inverse images 
under Hp of rank-one directions. 

Theorem 4.1 (Classification of Hamiltonians leading to elliptic Aronsson systems). 
Suppose that H g C^(IR^^") is a non-negative Hamiltonian with n > 1, N > 2. 
Consider the following statements: 

(i) There exists h G C^(§(M")),, h = h(p) with positive symmetric gradient hp = 
hp>0 such that 

(4.1) H{P) = h(^^P^P 

(ii) The Aronsson PDE system 

(4.2) Aoou := [Hp ® Hp + H [Hp]^Hpp^ {Du) : D^u = 
is quasilinear and degenerate elliptic, that is, the tensor map 

(4.3) ^»«(P) := Hp^^{P)Hp^^{P) + HiP)[Hp{P)]j^^Hp^^p^^{P) 
satisfies the strict Legendre-Hadamard condition and the symmetry condition 

(4.4) A{P) : {q ® w) ® {rj ® w) > 0, 

(4.5) A{P):{Q(g)R-R(g)Q) = 0, 
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for all TjeR^X {0}, w G M" \ {0} and P,Q,Re M^""" \ {0}. 

Then, (i) implies (ii). If moreover H is analytic at and satisfies 

(4.6) {Hp = 0} = {H ^0} = {0} , Hpp{0) > and Hpp : {v(»w-w(»v) = 0, 

forv,w G M", then, (ii) implies (i) when either 

a) n < 
or 

b) n > A and the q-th order derivative tensor Hp,,, p{{)) G (g)'''^-'(M"'^(K'M") is contained 
in the linear suhspace which consists of fully symmetric tensors T for which the 
only non-trivial components are of the form Ta-^ia2ja3k...agk, where am G {1, •■•,-^}; 
h3,k e {l,...,n}. 

If n > A but -ffp..._p(0) ^ , then H has the form (4.1) up to a forth order 
correction: H{P) = h{\P'^ P) + 0(1^1^^). 

In the case H{P) equals h(^^P^ P^ , the elliptic Aronsson system takes the form 

(4.7) AooU = (^Duhp <g) Duhp + h[Du\^<^hp^ : D^u = 

with h — h(^^Du^ Du^ . 

The extra assumption -ffp...p(0) G Ji"^ is necessary only in higher dimensions 
n > 4. It requires that Hp p{0) vanishes when more than 3 of its Latin indices are 
different to each other. The hnear space can be descrided as 

^« := {tg ®(«HM^®E") I T = r„,,,...„^,,ea,,, (8...®e„,,, : 

(4.8) . .Qz. ../3j. . . — . ./^-i. . .Qj?'.. . — ^...I3j...ai... ■} 

{ii,...iq} {i,j,k,...,k} =^ r = o|. 

If i?p...p(0) ^ ^''j then Hamiltonians with a httlc more comphcated fourth and 
higher order derivatives also give rise to elliptic systems. 

Proof of Theorem 4.1. We first prove the implication (j) (ii). For, assume that 
the Hamiltonian H has the form (4.1). We begin by observing that the symmetry 
assumption hp.. — hp., implies that second derivatives of h are fully symmetric in 
all indices: obviously since h is in C^(S(M")) we have hp..p^^ = f^pkipij ^^^'^ 

(4-9) ^PijPki ~ i^Pki)pij ~ i^pik)pij ~ ^PijPik- 

By using (4.9), we suppress the arguments in the notation of h and calculate 

(4.10) Hp^^P) ^ ^hp^^(Sc,f;Sk.,Pf}i + SafikiPfik) = Pakhp^^, 
and 

HPciPi3j{P) = 2^'''P^kPl„^{^M^ilP-frn + SpjSjmPjl^Pak + hp-^daf^Skj 

(4.11) = ^aphp^. + hp^^p.^PakPpm- 

Also, since hp > in §(M"), the operators Hp{Py , P^ : M" — > have the 
same nullspaces: 

(4.12) N{HpiP)'^) = {r/GM^ : i]'^ Php = 0} = N{P'^). 
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Hence, we obtain that [Hp{P)]^ = [P]^. By (4.10), (4.11), (4.12), we have 

{Hp®Hp + H[Hp]^HppYp) = Php®Php + h[P]^{l®hp + PhppP'^^^ 

(4.13) = Php ® Php + h[P]^ ^ hp, 

where h = /i(|P^P). Hence, in view of (4.3), equation (4.7) follows. Also, since 
h > and hp, [P]-*- are positive symmetric, conditions (4.4) and (4.5) follow as well: 

A{P) : {i] (g) w) (g) {rj (g> w) = Pakhp^^^rjaWiPpihp^^rjpWj + h[P]^priar]i3hp.^WiWj 

(4.14) = {Php : -q ® wY + h{[P]^ : ® vi){hp : w ® w) 

> 0, 

A(P) : iQ(g) R- R(g)Q) = ±Pakhp^,Qa^Pplhp,Ml3j ± /i[P]i^Q„»V.,^« 

(4.15) = ±{Php : Q){Php : R) ± h[P]-^ : (QhpR^) 

= 0, 

for aU 77 e M^, w e E", P,Q,R(E M^><". Hence, {ii) follows. Now we assume (ii) 
and prove the reverse implication. For, suppose H is analytic at and assume (4.4) 
- (4.8). By (4.5), we have 

(4.16) [Hp®Hp + H[Hp]^HppYp):{Q®R^R®Q) = 0, 

for all P,Q,Re M^^". By symmety of Hp® Hp and since by (4.6) we have H > 
and Hp ^0 on (M^><") \ {0}, (4.16) gives 

(4.17) [Hp]^Hpp:{Q®R-R®Q) = 0. 

By the identity [Hp]^ = I — [Hp]^ and since /, Hpp and [i?p]^ are symmetric, 
for Q = Ca® ei and R — ep ® Cj, (4.17) gives the commutativity relation 

(4.18) [Hp]i-,Hp^^p^^ = Hp^^pjHp]^^ 
on (E^><") \ {0}, that is 

(4.19) [Hp]^Hpp = Hpp[Hp]^. 

We set Aa,f;j := Hp^^p^^iQ). By assumtpion (4.6), we have A > in S(M^^"). By 
analyticity of H and since Hifi) = and -ffp(O) = 0, we have 

(4.20) H{P) = ^A:P®P + 0{\Pf), 

(4.21) Hp{P) = A:P + 0(|Pn, 

(4.22) Hpp{P) = A + 0{\P\), 

as |P| — ?> 0. Since A — Hpp{0) > and Hp{0) = 0, the map Hp is a diffeomorphism 
between open neighborhoods of zero in M^^". Hence, there is an r > such that 

(4.23) Hp : B^" := {Q e M^''" : \Q\ < r} Hp{M^") C 

is a diffeomorphism. Hence, there is a p > such that for < i < /?, ^ G §^^^ and 
w e §"~^, there exists a unique P{t) G B^" such that 

(4.24) ti®w = Hp{P{t)). 
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Moreover, \P{t) \ — > as t — > 0. The path P( . ) is the inverse image through Hp of 
the rank-one hne spanned by ^ (g) By (4.24), we have 

(4.25) [Hp{P{t))]'^ = [t((g}w]'^ = 

By evaluating (4.19) at P{t) and using (4.25) and (4.22), we obtain 

(4.26) (e®0(^ + o(l)) - (A + o(l))(e®0, 
at i ^ 0. In the Umit we get (^ (g) = A{(, (g) £,), that is 

(4.27) ^.u^kAkipj = AaiKj^K^/a, 

for aU i, j G {1, n}, a, /? G {1, N}. By the symmetry condition in assumption 
(4.6), for all i,j fixed the matrix Aaipj commutes with all l-dimensional projections 
of M^. Hence, it is simultaneously diagonalizable with them and as such a multiple 
of the identity. Thus, there is a symmetric matrix Aij such that 

(4.28) Aaifij — AijSap. 

Consequently, A : P(g>P = A : P^ P. We now set Baipj^k ■= Hp^^Pi^^p^,^{0). Then, 
by (4.28), equations (4.21) and (4.22) become 

(4.29) Hp{P) ^ pA + 0(|Pn, 

(4.30) Hpp{P) = 1(E) A +\b:P + 0{\P\^), 

and hence by (4.29) and (4.24) we get 

(4.31) t^(g>w = P{t)A + Oi\P{t)f). 

Since A > in S(]R^^"), we have i > in §(M") as weh. Thus, for < i < p, we 
have 

(4-32) ^ + = jj^^^^iiA-Yw). 

As t — >■ 0, we have |-P(i)| and by compactness along an infinitessimal sequence 
tm there exists a P with \P\ = 1 such that P{tm)/\P{tm) \ P- By passing to 
the limit in (4.32) as m — > cx) along {t,n}, we obtain that the limit of t,n/\P{tjn)\ 
exists and 

(4.33) Imr = P = 

m^oo \P(t,n)\ 



Since A-^ > and |P| = 1, for any v e §"-\ there is a u; G S""^ such that (4.33) 
becomes 

(4.34) hm = P = ^^v. 

By (4.19), (4.25), (4.29) (4.30), we have 



I®A+\B:P{t) + 0{\P{t)\^\ 



(4.35) = (l®A + \B:P{t) + 0(|P(t)n)c$5e 
By cancelling the commutative term ^ (g) ^(/ (g) A), (4.35) gives 

(4.36) + 0(^1))- (s:^ + 0(|P(t)|)) C ® 
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By passing to the limit in (4.3G) along t„i — > 0, in view of (4.34) we obtain 
(4.37) ^(E)^{B ■.£,(E)v) ^ {B ■.£,(E)v)^(E)^, 

for all C e S^-i, V e §"-1. Hence, (4.37) for v = Ck says 
(4.38) 

By (4.38), i? : ^ ® ^ is proportional to ^; hence, there is a map S : — > ®(^'E" 
such that B : ^ ® C -8(0 C, or 

(4.39) BatXjfj,kS,fj,^X = B^jkiOCa- 

By assumption (4.C) and induction, all second and higher order derivatives are fully 
symmetric in all their indices. Hence, we may fix fc € and suppress 

the dependence in them to obtain B^kX ^k^x = -B(C) with B e C°°(R^ \ {0}). 
The idea now is to differentiate in order to cancel both ^'s contracted with B and 
then contract again with a vector which annihilates ^ from the right hand side. For, 
by differentiating we get 

(4.40) DpBiO^^ = -B{£)5^p + 2B^p^^^. 

By (4.40), we obtain that DB{£^)(E)£^ is symmetric. Hence, we get that DB{^)(E)(^ — 
^ (g) Di3{^) and hence there exists B e C°°(M^ \ {0}), such that DB{£,) = B{£,)£,. 
Thus, (4.40) gives 

(4.4f) B(C)C®C +BiOl = 25 :e 

By differentiating the expression DB{(^) — B{£,)^, we get 

(4.42) DBiO(^^ = D^BiC) - B{^)I. 

By (4.42), we obtain that DB{£) ^s symmetric too. Hence, there exists B e 
C°°(M^ \ {0}), such that DB{1) = B{Cj£, and hence (4.42) gives 

(4.43) D^B{0 = Sm®^ + B{0I- 
By differentiating (4.40) again and inserting (4.43) we get 

(4.44) 2B^p^ = Dl^B{i)U + DpB{05^^ + D^B{S,)5^p 

for all ^ G M.^ \ {0}. Since N > 2, for each rj G we can choose a nonzero ^ 
normal to 77. Hence, by triple contraction in (4.44) we obtain 

(4.45) B:r,(Srj(8>fi = ^[Bi^i^v)^ + 5(01^1'] (^^^) = 0. 

Hence, by full symmetry in all indices we obtain Hp^iPf,jP^k{Q) — B^ipj^^ — and 
consequently third order derivatives vanish. We now set 

(4-46) Caipj^kSi := ^PciP/33P^fcP«(0) 
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and then for < i < p, (4.30) and (4.36) become 



(4.47) 
(4.48) 



Hpp{P{t)) ^I®A + -C ■.P{t)®P{t) + 0(|P(t)n, 



Pit) ^ Pit) 



omt)\) 



By setting t = tm and letting m oo, in view of (4.34), we get 



p(t)i \pm\ 



(4.49) 
for all ^ e 
(4.50) 



V e 



By (4.50), there exists a map C 
such that 



C : (C (g) w) (g) (^ «) u) 
I". Hence, for v = Ck, 

^\{0} 



C : (C fg w) fg (C ® ^') 



®(4)E" with Cy-fefc e C'°°(M^\{0}) 



(4.51) CaiKj\ktikS.K.^\^ti — CijkkiO^a- 

By fixing again the indices i,j,k, dropping them and arguing exactly as we did 
before for B^^^, there exist C,C e C°°{R^ \ {0}) such that 

(4.52) ilCa^p^s^s = CiOU^p^^ + CiOi^o^Sp^ + ^pS^^ + C^<5^„). 
By differentiating (4.52), we get 

(4.53) = C'(0(CaC,3'57<5 + £.P^jSaS + ^.j^aSps) 

Fix 77 G M^. Since TV > 2, there exists ^ ± ry, ^ 7^ 0. Then, (4.53) gives 

c(0 



3! 



(4.54) 

By (4.54), the function C is constant and moreover for all i,j,k, 

Cijkk 



VaVf3V-fV5 = 0{\r) ^1) 
= 0. 



(4.55) 



3! 



If either n < 3 or n > 4 but Hpppp{0) e where if* is given by (4.8), then 
in view of (4.46), the tensor Caipj^ykSi has no more than 3 different indices i,j,k,l 
for which it is non-zero. Hence, by full symmetry in all indices, (4.55) completely 
determines Hpppp{0) and gives 

(4.56) Hpppp{0) : (g(4)p ^ Ja,fc,P„,P„,P^feP/3, ^ 



(P^P) «) (P^P). 



Now we iterate the above arguments. The analog of (4.49) after blowing up along 
tm for q-th order derivatives is 



(4.57) ^®dHp...p{0) ■.®'-'^-^\^(giv) 



H 



P...P 



(0) 
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for all ^ e M^, V e M". When Hp,„p{0) e the only components of the tensor 
-ffpcjij ...Pa^i^ (0) which may not vanish are of the form 

(4.58) ifp.,.P„,,P„3....P»,.(0), 

where i,j,k e {l,...,n} and ai,...,aq £ {!,... ,iV}. Hence, (4.57), completely 
determines Hp p{0). By induction, all odd order derivatives of H vanish and all 
even order derivatives depend on P via P: we have 

J Cq : (^(9/2)pTp_ q ^ 2N, 
1 0, q e 2N+ 1, 



(4.59) Hp,„p{0) = 

q-th order 

for certain tensors C„ S (K-^'/^^M". Hence, by defining /i : M" (g) M" — ^ M by 



(4.60) Mp) 5] 2™C2™ : 

m— 1 

we obtain 

(4.61) H{P) - 

Hence, h>0 with /i e C°°(§(M")) and also hp = hj . Moreover by (4.3), (4.4) and 
(4.13), for w e §""1 and rj £ [P]^ D S^"\ we have 

(ARO\ h ( \ A(P) : (?? (g w) (g (?? (8 w) 

(4.62) h.p{p):w®w = , , , , , > 0, 

^ ' h{p)[P]-^ : T] (E) 7] 

where p = |/^^-P. Hence, hp > 0. If finally Hp,,p{0) ^ .if?, then iJ has the 
form (4.61) up to a correction of order OdPl**). This follows by decomposing each 
Hp,,,p(0) to the sum of a term in and a term in the orthogonal complement of 
.if' . The 0(|P|''') function arises from the series consisting of the forth and higher 
order parts of iJp...p(0) : (g)''^P in the orthogonal complements. The theorem 
follows. □ 

5. The 1- Dimensional Case: Aronsson's ODE System. 

5.1. Formal derivation of the general Aronsson ODE System. Let H he a 

non- negative Hamiltonian in C^(M x x M^), where N > 2 and we denote the 
arguments of H by H{x,ri, P). Consider the integral functional 



(5.1) Ejn{u,I) := J {H{x,u{x),u'{x))y"dx, 

where m > 2 and u : / C M — y M.^ . The Euler-Lagrange equation of functional 

(5.1) is the ODE system 

(5.2) (^H"'-\.,u,u')Hp{.,u,u'))' - H"''\.,u,u')H,{.,u,u') 
which after expansion and normalization gives 

(5.3) {H{_,u,u'))'Hp{_,u,u') + ^^^^^^[{Hp{_,u,u'))' -H,{_,u,u')) 

on / C M. We define the following projections of M^: 

(5.4) [i/(x,77,P)]T := sgn{Hpix,ij,P))®sgii{Hp{x,7],P)), 

(5.5) [H{x,r,,P)]^ := I - [ff (x, 77, P)]^. 



= 0, 
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Then, by employing (5.4) and (5.5) to expand the term in bracket of (5.3), we 
obtain 

(5.6) 

(ff (_, u, u'))'Hp{_, u, u') + ^^'"'"'^ [Hp{_, u, u')]^ (^{Hp{_, u, u'))' - iJ,(_, u')) 

- - "'^ [gp(-, u')]^({Hp{_, u, u'))' - g,(_, u')) . 

By perpendicularity of the orthogonal projections (5.4) and (5.5), the left and right 
hand sides of (5.6) are normal to each other. Hence, they both vanish. By re- 
normalizing the right hand side and rearranging, we get 

(i?(_, u, u'))'Hp{_, u, u') + H{_, u, u')[Hp{_, u, u')]^(^{Hp{_, u, u'))' - i7„(_, u, u')) 

(5.7) =- ^^^''_^'^\ hp{_,u,u')V ({Hp{_,u,u'))' - Hr,{^,u,u')). 

As TO — >■ oo, we obtain the complete system of Aronsson ODEs for a general Hamil- 
tonian with dependence on all the arguments 

{H{.,u,u'))'Hp{.,u,u') + H{_,u,u')- 

(5.8) ■[Hpi_,u,u')]^ (^{Hpi_,u,u')y - H.^i_,u,u')'^ ^ 0, 

whose solutions are curves u : / C M — > . 



5.2. Degenerate elliptic Aronsson ODE systems. We begin by observing that 
the Ellipticity Classification Theorem 4.1 readily extends to the case of H{x, rj, P) 
with dependence on all arguments; the form (4.1) of the Hamiltonian modifies to 

(5.9) H{x,rj,P) = /i(a:,7y,^F^p) 

and the PDE systems (4.2) and (4.7) modify by the appearance of first and lower 
order terns. In the case of ODEs where n = 1, the "geometric" Hamiltonians of the 
form (5.9) become the radially symmetric ones: 

(5.10) Hix,v,P) - h(^x,v,^\P\^y 

where h € C^(R x x [0,oo)) and the degenerate eUiptic Aronsson ODE sys- 
tem takes a particularly important and simple form. In the case of Aqo, we have 
h{x, r],p) = p. We now perform the derivation of the ODEs in the elliptic case. 

Suppose ft, G C^(M X x [0,oo)) with arguments denoted by h{x,ri,p) and 
define H e C^iM. x x M^) by means of (5.10). We henceforth assume 

(5.11) {hpix, r,, _) = 0} C {0} = {h{x, 77, _) = 0}, 

for all {x,ri) E M^+^. Assumption (5.11) is natural and will make the normal 
coefhcient H[Hp]^ of (5.8) continuous. By using (5.10) and supressing arguments, 
we calculate derivatives: 

(5.12) Hp = hpP, Hpp = hppP®P + hpl, Hr, = 

(5.13) Hpri = P^hprf, Hpx hpxP, Hx = h^. 
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By expanding derivatives in (5.8) and using (5.10), (5.12) and (5.13), we get 

{hpf{u'®u')u" + hp{u' ®hn)u' + h^hpu' 

(5.14) +h[hpu']-^(hpp{u' <Si u')u" + {u' ® hpr^)u' 

+ hpxu' + hpu" — = 0, 

where h = h[-,u, By assumption, (5.11), we have {hpu' = 0} = {u' = 0} = 

{h = 0}. Hence, we obtain that [hpu']-^ = [u']-^. On {u' ^ 0}, we multiply the 
normal term of (5.14) by to obtain 

{hpf{u'®u')u" + hp(^{u' ig) u')h^ + Ku'^ 

(5.15) + \u'\^hp[u']-^(hpu" - h,^ = 0. 
Hence, by using the identity \u'\^I — u' ®u' + \u'\^[u']-^, (5.15) gives 

(5.16) {hpf\u'\^u" - hp{\u'\^[l-2^^®^^h,, - Ku'^ = 0. 

By introducing the reflection operator : — > with respect to the hyper- 
plane [^]^, ^ e \ {0}, given by 




Figure 8. 

In view of (5.11), the systems (5.18) and (5.8) are equivalent on {u' = 0} as well. 
The system (5.18) comprises the degenerate elliptic Aronsson ODE system. 

Remark 5.1. We observe that in the special case where h — and h^ = 0, 

hx = 0, solutions of (5.18) trivialize to affine and actually (5.18) is equivalent 
to Aoo. In the special case where h — and ft.,, = 0, solutions of (5.18) 

become essentially scalar with affine rank- one range, that is u{M.) is contained in an 
affine line of since u" becomes proportional to u' and (5.18) becomes essentially 
scalar. Consequently, (5.18) is most interesting when h(^x,u{x), ||u'(x)p) depends 
on u(x) and hence h^j ^ 0. In this case the reflection operator R„' with respect 
to the normal hyperplane [u']^ is discontinuous on {u' = 0} at critical points of u, 
but the product ju'pR,;' is continuous. However, in all cases the system is always 
degenerate. 
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5.3. The initial value problem for elliptic Aronsson ODE systems. In this 
subsection we solve the initial value problem for ODE system (5.18) and consider 
some regularity questions. 

Theorem 5.2 (Tlie initial value problem for elliptic Aronsson ODEs). Suppose 
that h e C^(M X x [0,oo)) satisfies h,hp>0 and also (5.11) and consider the 
following problem for Aronsson ODEs 



AqqU = |u'p^/ipw" — R„'/i,j^ + h^u' — 0, 
u{xo) = Uo , u'{xo) = Wo , xo e M. 



(5.19) 
Then: 

(i) For any non-critical initial conditions (uoj'^o) G X (1^^\{0}), there exists a 
unique maximal smooth solution u : {xq — r, xo + r) — > M.^ for some r > which 
solves (5.19) and satisfies \u'\ > 0. 

(a) For any critical initial condition (mo,0) G M.^ x {0}, there exists at least one 
solution to (5.19), one of them being the constant one u = Uq. 

(iii) If 

(5.20) hri{x,ri,Q) ^ Q and hx{x,ri,p) — 0{p) as p 0, 

then bounded maximal solutions o/(5.19) starting (in positive time) from non- 
critical data, either are defined on [a;o,oo) being smooth and satisfying \u'\ > 0, 
or they reach a critical point u' = and form a discontinuity in u" in finite time. 

H If 

(5.21) c < hp < - forc>0, and hJx,rj,p) ^ Oip) as p — 0, 

c 

then bounded maximal solutions of (5.19) either are globally smooth or can be ex- 
tended past singularities as Wi^^{M.)^ strong solutions which satisfy (5.18) every- 
where and are eventually constant. 

The interpretation of W^^^ (M.)'^ solutions to (5.19) as strong everywhere so- 
lutions is the same as in Aronsson [Al, A2, A5]: at critical points of u whereon 
u" may not exist but is essentially bounded in a neighborhood of {u' = 0}, the 
coefficient vanishes. 

Example 5.3. The solution of problem (5.19) is generally non-unique for critical 
initial conditions. Choose h{x,rj,p) := +p. Then, (5.18) takes the form 



(5.22) I"T(w" - ^u'u) 







and the Hamiltonian is H{u,u') = 1(1^1 -f ju'l ). In view of example 3 in Aron- 
sson's paper [Al], for essentially scalar solutions u = £,v where ^ E E>^^^ and 
V : K — > M., (5.22) takes the form \v'\'^(v" -\- 1;)^ = 0. Hence, for initial conditions 
m( — I) = —ei, w'( — f ) =0, (5.23) admits the solutions ui{x) = eisina: and 
U2{x) = -ei. 

The non- uniqueness for critical data owes to that (5.18) is an 1-dimensional degen- 
erate elliptic system and the initial value problem is not well-posed for it. 

Proof of Theorem 5.2. All assertions follow directly by considering the following 
dynamical formulation of the ODE (5.18). For, we write the A^-dimensional second 
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order degenerate implicit system (5.18) as a 2iV-dimensional first order explicit 
system for a vector field defined off an A^-dimensional "slice" of R'^^ . For U = 

(5.23) U{x) := {u{x),u'{x))'^ , U : I QR — >R^^, 



(5.24) F{x,U) 



Ilvhjj{x,u, ^Iwp) 



/i^(£c,«,i|i;|^) 



:,«,i|«p) 

where 

(5.25) F : M X X (M^ \ {0}) — ^ R^^ . 

Then, in view of (5.23) and (5.24), ODE system (5.18) can be written as 

(5.26) tj{x) = F{x,U{x)) , U:ICR — >R^^. 
We now merely observe that the equation 

1 . / 1, „o\ hj_,u,^\u'\^) 



-u 



which follows by (5.26), implies that under assumption (5.20) the first term in the 
bracket becomes discontinuous at critical points of u, while the second one vanishes. 
Solutions extend past critical points where u" "jumps" by constant solutions. □ 
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